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The dynamical theory of thermally activated resonant magnetization tunneling in uniaxially 
anisotropic magnetic molecules such as Mni2Ac (S = 10) is developed. The observed slow dy- 
namics of the system is described by master equations for the populations of spin levels. The latter 
are obtained by the adiabatic elimination of fast degrees of freedom from the density matrix equation 
with the help of the perturbation theory developed earlier for the tunneling level splitting [D. A. 
Garanin, J. Phys. A, 24, L61 (1991)]. There exists a temperature range (thermally activated tun- 
neling) where the escape rate follows the Arrhenius law, but has a nonmonotonic dependence on the 
bias field due to tunneling at the top of the barrier. At lower temperatures this regime crosses over 
to the non- Arrhenius law (thermally assisted tunneling). The transition between the two regimes 
can be first or second order, depending on the transverse field, which can be tested in experiments. 
In both regimes the resonant maxima of the rate occur when spin levels in the two potential wells 
match at certain field values. In the thermally activated regime at low dissipation each resonance 
has a multitower self-similar structure with progressively narrowing peaks mounting on top of each 
other. [80163-1829(97)00141-0] 

PACS number(s): 75.45. -fj, 75.50.Tt 



I. INTRODUCTION 

In recent years there has been great experimental and 
theoretical effort to observe and interpret quantum tun- 
neling of magnetization in monodomain particles. The 
interest in this problem arises from the fact that the 
magnetization M of a particle containing a few thousand 
atoms is a macroscopic degree of freedom. Thus tunnel- 
ing of the particle's magnetization between different equi- 
librium orientations at low temperatures requires strong 
coherence between atomic spins and may be very sensi- 
tive to the interaction with the environment. A similar 
problem has been extensively studied in superconductors 
in the context of macroscopic quantum tunneling, where 
good agreeiaent has been achieved between theoryld and 
experiments Observation of magnetization tunneling is 
complicated by the difficulty in preparing identical wag- 
netic particles. Experiments have been performedu on 
particles distributed over sizes and shapes. These exper- 
iments revealed temperature-independent magnetic re- 
laxation which was attributed to tunneling. When an 
effort was rnadc to narrow the distribution, resonance 
was observedalj in the absorption of the ac field, similar 
to the tunneling resonance in the ammonia molecule. 

Difficulties in manufacturing identical magnetic par- 
ticles for tunneling experiments have led Xa, new 
techniques of measuring individual particleaS'EI and 
to the idea of searching for magnetization tunneling 
in magnetic molecules of large spin. The system 



that caught the most recent attention is the crystal 
Mni2 acetate (Mni2Ac) having the chemical formula 
[Mni20i2(CH3COO)i6(H20)4]-2CH3CQOH-4H20. This 
compound has been synthesized by Lisp but its physi- 
cal properties had not received much attention until Ses- 
soli et aln noticed magnetic bistability of this system. 
In the Mni2Ac molecule the 12 Mn ions are strongly 
bound ferrimagnetically via the superexchange through 
oxygen bridges. These molecules behave effectively as 
magnetic clusters of spin S = 10,0 as has been confirmed 
by the Curie-law temperature dependence of the suscep- 
tibility X- As follows from the very low value of the 
Curie constant, Oc ~ —0.05 KjlJthe interaction between 
the Mni2Ac molecules is very weak, presumably of the 
dipole-dipole origin. Mni2Ac is characterized by a very 
strong uniaxial anisotrapy TYa = ^DS^, where D ~ 0.72 
K from high-field EPR,L1 D ~ 0.75 K from single-crystal 
magnetic susceptibilitytj measurements—and D ~ 0.77 
K from neutron scattering experiments. Ilj This leads to 
a barrier of about U = DS^ ~ 75 K between the states 
±S. Notpy, however, that experiments on resonant spin 
tunnelingEJ (see below) suggest a value of D close to 0.6 
K and correspondingly the barrier height of 60 K. 

The advantage of Mni2Ac and other molecular mag- 
nets is that they are rather simple model systems, which 
facilitates their theoretical consideration and interpreta- 
tion of experiments. Of course, it should be understood 
that a cluster of spin 10 cannot be treated macroscopi- 
cally. The limit of macroscopic quantum tunneling is the 
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one where the quantization of spin levels is irrelevant. 
On the contrary, in the Mni2Ac cluster the distance be- 
tween the ground-state and the first excited level is 12-15 
K. At low temperature quantization of levels must, there- 
fore, dominate the properties of the system. In this sense 
Mni2Ac is closer to conventional quantum-mechanical 
systems where tunneling is of a resonant character. Nev- 
ertheless, as we shall see, the high value of spin leads 
to the macroscopic time scale for the dynamics of the 
magnetization, which has been tested in macroscopic ex- 
periments. 

An important feature of Mni2Ac is that if no strong 
transverse field is applied to the system, the interactions 
responsible for tunneling are small in comparison to the 
anisotropy energy TYa = —DS^ which itself conserves the 
Sz component of the spin. As a result of this and of the 
large spin of the system, the tunneling between low-lying 
energy levels should be extraordinary slow, which makes 
Mni2Ac an excellent candidate for information storage 
at the molecular level. Another possible application of 
molecular magnets is that for quantum computing. For 
that application tunneling between the low-lying states 
should be made more pronounced, and the interaction 
with the environment destroying coherent oscillations of 
the spin between two wells should be kept small. This 
is hardly the case for Mni2Ac where nuclear spias of 
manganese atoms strongly suppress the coherence. Il3 The 
example of Mni2Ac is, however, instructive since other 
systems with similar properties can be developed, which 
could be better candidates for quantum computation. 

The first indications of magnetization tunneling in 
Mni2Ac were seen in the magpetization relaxation ex- 
periments of Paulsen and Parkllj and the dynaroii; sus- 
ceptibility measurements of Novak and Sessoli.Ej The 
measured relaxation rate of Mni2Ac followed the Arrhe- 
nius law T — Fq exp(— f/es/T) with the peaks at some 
values of tlie longitudinal field Hz- These peaks were 
interpretedll3 as the resonant thermally assisted tunnel- 
ing between the levels near the top of the barrier , which 
decreased the effective barrier lieight UcS. Subsequent 
dynamic hysteresis experimentsEj have proved that con- 
jecture as they have shown many regularly spaced steps 
in the hysteresis loop at the values of at which the 
levels on both sides of the barrier come into resonance 
(see also Refs. |l^, |l8|, |^ |o|). These steps indicate an 
increased relaxation rate at the corresponding bias fields 
Hz- Very recently a similar observation was made on 
Mni2 phosphatO which was described as a magnetic clus- 
ter of spin S — 9.5. 

The transverse field H^ applied to a uniaxial mag- 
netic system mixes the unperturbed energy levels and 
enhances tunneling. The search for an increased tunnel- 
ing in thcj-transverse field has been undertaken in recent 
hysteresisE3 and dynamic susceptibilityBj measurements. 
The results show that the speeding up of the relaxation 
can be explained mostly through the classical effect of the 
barrier lowering in a transverse field, whereas the reso- 
nant tunneling peaks remaining after subtraction of this 



main effect are nearly independent of H^- Actually both 
effects come from the same source: The classical height 
of the barrier can be determined quantum mechanically 
from the condition that the tunneling level splitting be- 
comes comparable with the level spacing, which means 
strong nonresonnat tunneling, i.e., the absense of a bar- 
rier at that level Eil 

A large number of experimental observations of mag- 
netization tunneling in molecular magnets has been ac- 
cumulated to date and the major relevant physical pro- 
cesses have been identified. A theoretical framework for 
the dynamical description of the combined process of the 
thermal activationand tunneling in these materials is still 
lacking, however. In particular, the form and the width 
of the tunneling peaks measured in experiments has not 
yet been explained. The aim of this article is to supply 
an appropriate theory. 

The idea of the work is to apply the density matrix 
formalism in the case when the tunneling is caused by a 
transverse field H^ which is small enough and can be con- 
sidered as a perturbation. The applicability criterium of 
this method is H^ <C Hp^, where Hp^ = {2S — l)D/{g^B) 
is the anisotropy field. The latter in turn coincides with 
the critical value of the transverse field at which in the 
classical case of 5 ^ 1 the double-well structure of the 
spin energy disappears. For Mni2Ac, the anisotropy field 
is of order 10 T, so that the condition H^ ^ Ha allows 
for rather large H^ - In this relevant range of the trans- 
verse field one can use the physically transparent and 
technically convenient basis of the eigenfunctions of the 
anisotropy energy Ha = —DSl- The slow dynamics of 
the system driven by the thermal activation and tunnel- 
ing processes can be described with the help of the adi- 
abatic elimination of the fast degrees of freedom in the 
density matrix. The latter is a dynamical generalization 
of the calculation of the tunneling leveLsplittings in the 
high orders of the perturbation theory.E3 

The remaining part of the paper is organized as fol- 
lows. In Sec. |lj the properties of an isolated magnetic 
cluster in a transverse field are briefly reviewed and the 
perturbation theory is compared with other approaches 
to the problem. In Sec. Ill the density matrix equation 



(DME) for the uniaxial magnetic system interacting with 



a phonon bath is formulated and discussed. In Sec. IV 
the fast degrees of freedom in the DME are eliminated 
and a simplified system of equations describing the slow 
spin dynamics in terms of the diagonal and antidiago- 
nal matrix elements connecting resonant pairs of levels 
in different wells is derived. It is shown that the level 
broadening due to the interaction with the environment 
suppresses coherent oscillations and, if strong enough, 
makes the motion of the spin between two degenerate 
levels overdamped. In this case, and also in the case 
of the thermally assisted quantum tunneling, when the 
relaxation rate is limited by the exponentially slow pro- 
cess of climbing up the energy barrier, the DME further 
simplifies to the system of kinetic balance equations for 
the level populations iV„i only. The latter describes the 
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hopping of particles between adjaicent energy levels and 
through the barrier. In Sec. the system of equations 
for the level populations Nm is solved analytically in the 
Arrhenius regime T <^ U DS^ . In Sec. VI the tran- 



sition from the Arrhenius regime to pure quantum tun 
neling at lower temperatures is discussed. In Sec. VII 



the numerical results for the dependences of the escape 
rate on longitudinal and transverse fields in the Arrhenius 
regime are presented. Here we also analyze the influence 
of the Mn nuclear spins and a small scatter of the easy- 
axis directions in the oriented pol ycrys talls on resonant 
magnetization tunneling. In Sec. VIII further develop- 
ments of the theory and suggestions for experiments are 
discussed. 



II. TUNNELING LEVEL SPLITTING AND 
CLASSICAL BARRIER LOWERING 

The spin Hamiltonian of an isolated Mni2Ac molecule 
in magnetic field H can be written in the form 
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-DS, — HzSz — HxSx 



(2.1) 



where H stands for g/isH with g ~ 1.9. Henceforth we 
will usually drop the combination gfj,B for better read- 
ability of the formulas. The system is described by the 
25' -I- 1 energy levels which in the absense of the trans- 
verse field are labeled by the spin projection m on 
the z axis and given by Sm — —Dm? — Hzm (see Fig. 
|l|) . It can be easily checked that for the regularly spaced 
values of the longitudinal field Hz satisfying 



Hz^Hzk^kD, fc = 0,±l,±2,. 



(2.2) 



the energy levels on both sides of the barrier are pairwise 
degenerate 



TO < 0, 



(2.3) 



The latest high-field EPR experimentsHj suggest that 
there are correction terms of the types —ASj and 
-B{Sl + Si) in the spin Hamiltonian (|1|) of MuizAc. 
This means that the degeneracy of different level pairs 
TO, to' is actually achieved at slightly different values of 
Hz- We shall, however, ignore this effect in the following 
since it does not significantly change the results. As we 
shall see, only one or maximally two pairs of degenerate 
levels contribute to resonant tunneling, and hence the 
lack of simultaneous degeneracy of all appropriate level 
pairs is unimportant. 

The model Hamiltonian ( |2.lD was a whetstone for dif- 
ferent theories of spin tunneling long before its relevance 
for Mni2Ac and other molecular magnets had been es- 
tablished. In the quasiclassical limit S ^ 1, the rate 
of tunneling from the ground-state for different values of 
Hz was calculated by Chudnovsky and GuntheiO with 
an exponential accuracy with- the help of the instanton 
technique. Enz and Schillingc3 developed a more sophis- 
ticated version of the instanton approach to spins to ob- 
tain the ground-state tunneling level splitting with the 
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FIG. 1. Spin energy levels of a Mni2Ac molecule for Hx = 
and Hz~ D corresponding to the first resonance, fc = 1, in 
Eq. 01). 



prefactor. The latter result was rederived by Zaslavskiica 
by a more simple method based on the mapping, pato a 
particle problem. Also, van Hemmen and SiitcE^'Ell for- 
mulated the WKB method for spin systems and calcu- 
lated the tunneling rates and cor resp onding level split- 
tings for the excited states of Eq. ( ^.l|) . Scharf, Wreszin- 
ski, and van Hcmmenc3 proposed an approach based on 
a particle mapping with subsequent application of the 
WKB approximation to refine the results for the split- 
tings of excited levels for systems with moderate spin. 
The applicability of this approach is confined, however, 
to the limit of small transverse fields H^ , where it is still 
possible to label the energy levels of Eq. (2.1) by the 
quantum number m. 

In the case of small H^, which, as we shall see, is rele- 
vant for magnetic clusters with moderate spin, the level 
splittings can be calculated in a more direct and sim- 
ple way using the high-order perturbation theory. An 
early application of this method is due to Korcnblit and 
ShendeiH who studied ground-state splitting in rare- 
earth compounds having high spin values (e.g., S = 8 for 
Ho). GaraninLJ has derived a formula for the splitting 
of all levels of the Hamiltonian (2.1). A recent remake 
of the method is due toJIartmann-Boutron.O Schatzer, 
Breymann, and ThomasEs extended the perturbative ap- 
proach to describe tunneling in a system of two spins. 

In the general biased case, the tunneling level split- 
ting of the resonant level pair m, m' appears, minimally. 
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FIG. 2. Tunneling splittings Asmrn' for Hz — and dif- 
ferent values of the transversa field. The results of Scharf, 
Wreszinski, and van Hemmeno are indiscernible from the per- 
turbative ones in this scale and they are shown only for the 
one value of the transverse field. 



in the |m — m'|th order of a perturbation theory and is 
given by the shortest chain of matrix elements and energy 
denominators connecting the states m and m', 

xKn+l,m+2 

where 

Ki,m+1 = {'m\H,j.Sx\m + 1) = -Hxlm,m+li 



1 



£m+2 — £71 



Vm'-l.' 



(2.4) 



lTn,m+i = y/ S{S + 1) — m{m + 1) are the matrix ele- 
ments of the operator Sx, which are symmetric functions 
of their arguments, and Em = -^Dm? — Hzm are the un- 
perturbed energy levels. The calculation in Eq. (2.4VJpr 



the arbitrary resonance number k yields the formula 
2D 



[(to' — to — 1)!]^ 
' {S + m')l{S - my. (H, 



{S - m')\{S + m)\ \2D 



(2.6) 



which is the generalization of the zero-bias result of Ref. 
P5|. Note that here, according to the convention of Eq. 
(pjl) , TO < 0, to' = —m — k, and hence m! ^ra > 0. Equa- 
tion (2^) describes the interaction between the pair of 
resonant levels m, to' through the intermediate levels in 
the virtual state. As is well known for the two-state prob- 
lem, the splitting Aemm' is exactly equal to the tunneling 
frequency f2mm' with which the probability of finding the 
system in one of these states oscillates with time if the 
initial condition is an unperturbed eigenstate. 



The tunneling splittings given by Eq. (^^) are repre- 
sented in Fig. g for Hz — and different values of the 
transverse field, in comparison with the results of other 
approaches. One can see that the splittings change by 
orders of magnitude with changing to by 1. If the split- 
ting of the pair to, to' becomes comparable to the level 
spacing in the well, which is of order 2Z3|to|, the tunnel- 
ing becomes strong and of nonresonant character; i.e., the 
barrier for a particle going into the other well disappears. 
For this pair of levels the perturbation theory clearly 
breaks down, but for the next lower pair to — 1, to' 4- 1 
(see Fig. |^) it already works well. 

The sharp boundary between the levels localized in 
one of the wells and the delocalized ones, which was ob- 
served above, is also characteristic for the classical theory 
where there is a similar separation between the localized 
and escape orbits ai-some energy. Accordingly, as was 
shown by Friedmanj23 the transverse-field dependence of 
the classical barrier height. 



Hx 
2SD' 



(2.7) 



can be reproduced fo r sm all hx with the help of the 
perturbative formula (2_^). Indeed, in the quasiclassi- 
cal hmit 1 < |to| < 5 Eq. for = can be 

simplified to 



ASr. 



2D\m\ fHxSe' 



2\m\ 



(2.8) 



and compared to the level spacing D\m\ to obtain the 
value TO = TTif) at which the barrier is effectively cut by 
the tunneling. For \m\ ^ 1, the value of TOb can be found 
with a good accuracy by equating the fraction in brackets 
to unity. The result has the form 



(2.5) inEq. (^ 



TO? = 2S^hx — , 



(2.9) 



which leads to the effective barrier height U = DS^ — 
Drni = DS'^[1 — 2/ia;(e^/8)]. This is in accordance with 



Eq. ( ^?7| ) for < 1, except for the factor e^/8 « 0.92. 
The nontrivial feature of this derivation is that the re- 
sulting classical barrier lowering is of first order in hx^ 
although the corrections to the energy levels arise only 
in the second order of the perturbation theory. The latter 
have the form 



.(2) 



-Dm'' I 1 



2S'2[S'(5+1) 



^2 [4^2 _ j^] 



(2.10) 



It can be checked that for m = TOf, and 1 <C mi, <^ S 
the correction term in the curly brackets makes up the 
universal number 8/e'* « 0.15. This means that near 
the renormalized top of the barrier , to = to^, the per- 
turbation theory relies on a small numerical par ameter 
rather than on hx- The artifact e^/8 in Eq. (2.9) is the 
consequence of dropping the effect of the level mixi ng in - 
side the wells described to lowest order by Eq. (2.10). 
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It should he .noted that e^/8 also appears in the WKB 
resultfCjOH for the tunneling level splitting in the case 
of small transverse field, and it can be attributed to the 
inaccuracy of the WKB method near the top of the bar- 
rier. In this article, we will neglect these effects and 
study the tunneling transitions between the wells per- 
turbatively in the basis of the eigenstates of the operator 
Sz- It should be noted in addition that, as was checked 
by Chudnovsky and Friedman,c3 the level matching con- 
dition ( ^.2| ) remains uneffected by the transverse field at 
least up to fourth order in h^- 

Now let us consider the question how the level splitting 
changes from one level pair to another in more detail. For 
the pairs of resonant lev els shown in Fig. |l|, with the use 
of the basic formula (2_^) in the unbiased case, one comes 
to the result 



Ae. 



m+l,m' — 1 



2|L|) |m| 



(2.11) 



where rrih is given by Eq. (2.9). One striking implica- 
tion of this formula is that the splitting ratio is large 
everywhere in the wells: Even near the top of the renor- 
malized barrier, m ~ to;,, the tunneling splitting changes 
by a large factor e"* « 55, moving one step up the barrier. 
This universal behavior, independent of the spin value S 
for S* ^ 1, shows that even in the quasiclassical limit 
the tunneling splitting cannot be treated as a smooth 
function of the energy. The determination of the level at 
which the barrier disappears is, th erefore, quite precise. 
Another consequence of Eq. (2.11) is that resonant tun- 
neling is to the same extent inherent in models of large 
spin S as in those of moderate spin. 



III. SPIN-BATH INTERACTIONS AND THE 
DENSITY MATRIX EQUATION 

The thermally activated escape of the Mni2Ac spin 
over the potential barrier DS^ ~ 70 K is accompanied 
by the transitions between the energy levels with the en- 
ergy differences ranging from D{2S — 1) ~ 13 K near 
the bottoms of the potential wells to D ~ 0.7 K near 
the top of the barrier. Such a process requires an energy 
exchange between S and other degrees of freedom of the 
whole system. 

The dipole-dipolc interactions between different mag- 
netic clusters contribute to the macroscopic magnetic in- 
duction B = H + 4ttM which is actually "felt" by the 
spins and which should replace the external field H in all 
the formulas for spin tunneling and thermal activa ti o n ,. 
As was shown in dynamic hysteresis experiments Jia't^l 
this internal field correction is quite essential for a careful 
analysis of the experimental data. The fluctuating part 
of the dipole-dipolc interactions which could cause the 



spin relaxation iias been shown to be inefficient by dilut- 
ing the sample.E3 Indeed, this interaction is of the order 
of the dipole-dipole energy of two neighboring clusters, 
— [gt^sSY /vq, where vq is volume of ihe unit cell. 
Using g = 1.9 and = (17.31)^ x 12.4A,eI one obtains 
Ed — 0.06 K [in accordance with the measured value of 
the Curie constant Qc = -0.05 K (Ref. ^] which is 
much smaller than the distances between tlie energy lev- 
els. There is also a more subtle argument ,eZI based upon 
energy conservation and the nonequidistant character of 
the spin energy levels, which rules out the contribution 
of dipole-dipole interactions to the relaxation in the tem- 
perature range T <C [/. 

The nuclear subsystem also cannot supply energies 
which would be large enough for the relaxation over the 
70 K barrier. Nevertheless, nuclear spins produce a hy- 
perfine field on the effective electronic spin, which can 
give rise to tunn eling . This mechanism will be considered 
in detail in Sec. VII. Here we will describe tunneling as 



caused by the externally applied transverse field Hx- 

The remaining two types of the interaction of a 
Mni2Ac spin with the environment are those with 
phonons and photons. Unlike the interactions reviewed 
above, the phonon and photon subsystems play the role 
of a thermal bath, rendering the spin subsystem a defi- 
nite externally controlled temperature. It can be imme- 
diately seen that in the presense of phonons the photon 
processes can be safely neglected, since the light velocity 
c is much greater than the sound velocity v and, as a 
result, the photon density of states is smaller than the 
phonon one. At low temperatures the leading processes 
are the emission and absorption of phonons, accompa- 
nied by the hopping of spin between energy levels. At 
higher temperatures Raman scattering processes can be- 
come dominant. The energies of phonons in Mni2Ac are 
large enough for the exchange with the spin-subsystem: 
As follows from specific heat measurementsjlj the Debye 
temperature 9d corresponding to the phonon energy at 
the edge of the Brillouin zone is about 36 K. 

Spin-phonon interactions in materials with a strong 
crystal-field anisotropy are mainly due to the modula- 
tion of the crystal field by phonons—This mechanism was 
extensively studied in past years .c3 The possible spin- 
phonon coupling terms for substances of different symme- 
tries are listed in Ref. For Mni2Ac and other molec- 
ular magnets, the spin-phonon interactions, as well as 
the (presumably complicated) phonon modes themselves, 
have not yet been investigated. Moreover, an attempt to 
describe the interaction with phonons rigorously would 
lead to a serious complication of the formalism without 
bringing any new qualitative results. We will resort to 
various simplifications, assuming, in particular, that the 
phonon spectra of molecular magnets contain, as for an 
isotropic elastic body, one longitudinal and two trans- 
verse modes. Similar simplifications were also made in 



Ref. 37, where the pure thermal activation escape rate in 
Mni2Ac was studied. 

The lowest-order spin-phonon interactions allowed by 
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the time-reversal symmetry are linear in phonon opera- 
tors and bilinear in the spin operator components, con- 
taining various combinations S^S^, where a, (3 = ±,z. 
The simplest of these interactions is due to the, rotation 
of the anisotropy axis by transverse phonons.E3 We will 
use this mechanism for the illustration of our method 
since it does not employ any unknown characteristics of 
the crystal-field distortions accompanying other types of 
lattice vibrations. 

For the arbitrarily oriented anisotropy axis n, the 
anisotropy part of the spin Hamiltonian (2.1) can be writ- 
ten as Ha = —D{nS)^. Transverse phonons change the 
vector n by (5n = [54> x n], where Scf) — (1/2)V x u is 
the local rotation of the lattice and u is the lattice dis- 
placement. The first order term on Sn in Ha gives the 
spin-phonon Hamiltonian which in coordinate form reads 



Hsp = D{Sz,Sx}uJzx + D{Sz,Sy}ujz 



where 



2 \drR dva 



(3.1) 



(3.2) 



and {Sa, Sp\ is the anticommutator. In terms of phonon 
operators ak\ and a\y^^ 

where M is the unit cell mass, N is the number of cells 
in the lattice, CkA is the phonon polarization vector, A = 
t, t, I is the polarization, and iOkx = v\k is th e ph onon 
frequency. Performing differentiation in Eq. ( |3.2| ), one 
can transform Eq. (|3.lD to 
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Ari/2 ^ 

kA 



V^fe{5„ (r,kAS)}(afeA 



"kXl 



Here the spin-phonon amplitide Vk is given by 



D ( U!kt 



23/2 I a 



1/2 



rif = Mvf 



(3.4) 



(3.5) 



and the vector rjkA is determined by 



the thermoactivation process in uniaxial spin systems, as 
Mni2Ac, in Refs. ^ and If a transverse field or an- 
other level mixing perturbation is applied to the system, 
the nondiagonal elements of the DME appear, whose slow 
dynamics describes the tunneling process. The major ad- 
vantage of the DME is that it provides a natural account 
of resonant tunneling in systems of moderate spin, which 
is lost in quasiclassical approaches for truly macroscopic 
systems. 

A common routine for obtaining a system of kinetic 
balance equations is to calculate the transition proba- 
bilities according the Fermi golden rule and then insert 
them into the equations that are themselves postulated 
but not derived. Such an approach is methodically insuf- 
ficient since the transition probabilities are obtained with 
the help of the time dependent perturbation theory where 
the probability of finding the system in states differing 
from the initial fully occupied state are used as a small 
parameter. In other words, this method describes only 
the initial stage of the relaxation process for a special 
type of initial conditions. Although it incidently leads to 
the correct master equation, the same is not true for the 
general DME. Indeed, spin-phonon couplings of the type 
J2k ^fc'S'f flfeO^, corresponding to the elastic scattering of 
phonons, do not result in transitions between the energy 
levels and do not contribute to the coefficients of the mas- 
ter equation. On the other hand, such terms modulate 
the energy levels and contribute to the linewidths, which 
manifest themselves in the dynamics of the nondiagonal 
elements of the density matrix. 

A rigorous method of the derivation of the density ma- 
trix equation valid fm all.iiines employs the projection 
operator technique .l£jtijcjl£j For spin systems, the de- 
tails of calculations are described in Ref. The re- 
sulting DME can be found in Ref. where the model 
without single-site anisotropy, accounting for both one- 
phonon and Raman scattering processes, was used to 
derive the Landau-Lifshitz-Bloch equation for ferromag- 
nets. This DME is written in terms of the Hubbard oper- 
ators X™" = \m){n\ forming the complete basis for the 
spin subsystem. In the Heisenberg representation the 
operators X™" are related to the spin density matrix: 
p„in = (A r'""(t )). For the present model described by 
Eqs. (|2]l|) and the resuhing DME reads 



= 0, 



x.y 
»/kA 



z x.y 
SkA"k 



ekA"i< 



(3.6) 



where nk = k/fc. On can see that the coupling to lon- 
gitudinal phonons in Eq. (3.4) vanishes, as it should be, 
since eki = nk. 

The evolution of a spin system coupled to an equilib- 
rium heat bath can be described by the density matrix 
equation. The diagonal elements of the density matrix, 
Pmm = Nm, describc the population of the energy levels. 
In the absense of interactions noncommuting with Sz in 
the spin Hamiltonian Ti, the DME reduces to the closed 
system of kinetic balance equations, or master equations, 
for the populations Nm in the basis of the eigenstates 
of the operator Sz- The latter was applied to describe 



\rmn . \rmn rj ( 1 -v^m+l,n i ; x^m— l.n 



m,n+l 



(3.7) 



where uJmn = £m — £n are the frequencies associated with 
the transition n —> m, the unperturbed energy levels e„ 
are given by Em = —Dw? — HzUi, the factors gjiB and Ti 
are dropped for convenience, the matrix elements lm,m±i 
are given by Eq. (2.5), and Rmn is the relaxation term. 
The latter has the non-Markovian form 



fdt'^Y.^,^Afkit'-t)-Bfkit 



(3.8) 
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where 



A = Qs{t')[Qsit),X^^{t)], 
B=[Qs{t),X"^^\t)]Qsit'), 



(3.9) 



the spin operator combination Qs = {Sz, ('7kAS)} conies 
from the spin-phonon Hamiltonian ( |3.4| ), the function 
fk (t) characterizing the bath in the present case of the 
one-phonon processes is given by 



(3.10) 



Uk = (e"*"/-^ — 1)^^ are the boson occupation numbers, 
and uJk = t okt a rc the frequencies of transverse phonons. 

In Eq. ( |3.9| ) the spin operators should be expanded 
over the X™" basis as foUows: 



s-i 



S± = S,±lSy. (3.11) 



s. 



m=-S 
S-1 

ra— — S 



For th e o ne-phonon processes, the integral over t' in 
Eq. (B.8) converges on the scale of l/ujmn which is 
much shorter than the relaxation time of the spin system. 
Hence, the lower limit of this integral can be extended to 
to = — oo and the t' dependences of the operators X'"" in 
the relaxation term can be considere d as governed solely 
by the conservative part of the DME (3/7). Finding these 
time dependences is a matter of numerical work, if the 
transverse field Hx is not small. Here serious complica- 
tions arise, since the evolution of each operator X™" is 
a linear combination of all possible types of spin motion. 
This means simply that the unperturbed basis we have 
chosen is not suitable in situations with strong level mix- 
ing. However, in the case of small Hx one can neglect 
these effects and use the unperturbed time dependences 



(3.12) 



N ow one can calculate combinations A and B in Eq. 
( |3.9|) with the use of the representations ( 3.11 ) and the 



equal-time relation x™" X = X'^'^Ski which replaces 
the commutation relations for the spin components. The 
sum over the ph ono n polarizations A in Eq. (3.8|) can be 



done using Eq. (3.6) and the property of the polarization 



vectors X^a ^a^a — ^ap- Neglecting the imaginary part 
of the relaxation term Rmm corresponding to the renor- 
malization of the spin energy levels due to the coupling 
to the bath, one arrives at the final form of Rmn- 

R + WnM+l]X"^^ 



+ ■^lm,m-lln,n~l[Wm,m-l + Wn,n-l\X"'' 



Here 1^ ,^±i = Zm,m±i(2TO ± 1) with the factor 2m ± 1 
coming from the operator Sz in Eq. (3^), and the uni- 
versal rate constant Wmn = W{uJmn) of the one-phonon 
processes is given by 



W{uj) 



V^{{nk + l)TrS{uJk+i^) 
'nkTr6{ujk-u!)}, (3.14) 



where vq is the unit cell volume and the overall factor 
2/3 says that only two transverse modes of the total 
three phonon modes are active in the relaxation mech- 
anism under consideration. One can check that the rate 
constant satisfies the detailed balance condition W{lu) — 
W{—uj)exp(—uj/T). At low temperatures phonons die 
out and W{uj) with w > 0, which corresponds to the ab- 
sorption of a phonon, becomes exponentially small. The 
result for W{uj) with uj < (the emission o f a phonon) 
calculated with the help of Eqs. ( 3.14 ) and (3^) reads 



Here we have used 9^ 



(cf. Ref. |9[ 

the Debye temperature 9o ~ fiiOk 



{hvtf/vo for 
The constant 6 is 



defined as = ^t^D — ^"^P^^fi where p is the density 



and Qt is given by Eq. (3.5). 

Note that Eq. (O with R^n given by Eq. (|3.13|) 



is still an operator equation, and the equation of motion 
for the density matrix elements, p,„„ = (X™"), should be 
obtained by taking its quantum-statistical average over 
the initial state of the spin. This is, however, a trivial 

task, since the equation for X"*" is linear. 

In the c ase = the density matrix equation (3.7) 
and (3.13) reduces to a system of kinetic balance equa- 
tions for the diagonal elements Nm = X"™ , the equilib- 
rium solution of which is given by 



1 



s 

E 



-e™/T 



(3.16) 



The thermoactivation relaxation rate F in the model with 
Hx — was studied in Ref. ^ and recently in Ref. 

In the latter work Raman scattering processes have 
also been taken into account, and the spin relaxation 
rate was calculated for arbitrary ratios U/T in terms of 
the integral relaxation time Tint. It was shown that in 
systems with larger spin values, even in the Arrhenius 
regime U/T ^ 1, there are several limiting cases for the 
prefactor Fq in the expression F = Fg exp(— f7/T) as a 
result of the interplay between the one-phonon and Ra- 
man scattering processes. Here we concentrate on the 
low-temperature region, and thus only one-phonon pro- 
cesses will be considered. 



+ ?^.„_iVF„_i,„]X"". (3.13) 
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IV. SLOW DYNAMICS OF THE DENSITY 
MATRIX: COHERENCE AND TUNNELING 
BETWEEN RESONANT LEVELS 

The possible frequencies, with which the density ma- 
trix elements X™" evolve in time according to the DME 
(3.7), range from ujqs — DS^ (for = 0) to very small 



ones corresponding to overbarrier relaxation and tunnel- 
ing. In the low-temperature range T ^ U, these fast 
motions decay with the rate corresponding to the relax- 
ation inside one well, which is much larger than the ther- 
moactivation escape rate or the tunneling rates. In the 
long-time or low-frequency dynamics, the variables X™" 
corresponding to the large ujmn play the role of "slave" 
degrees of freedom, adjusting themselves to the evolution 
of the slow variables, and hence they can be adiabatically 
eliminated. 

The slow variables of our problem are the diagonal 
matrix elements Nm = X™™" , as well as the antidiagonal 
elements X"™ whose transition frequency ujmm' = £m ~ 
Sm' is the detuning of the resonant levels m and m' 



= {H - Hk){m' - m) 



(4.1) 



[cf. Eqs. (2.2) and ( |2.3| )]. The equations of motion for 
these slo w va riables can be obtained in the following way. 
In Eq. (1^) for X""™, the terms containing 
and X™''"''""'", which are generated minimally by nonzero 
X™ ™ and X™™ , correspondingly, are responsible for 
tunneling in the lowest approximation. In the dynamical 
equations for these elements one can neglect the terms 
j^m+i,m g^j^j j^m,m.+i^ .^gjj l-]^g relaxation terms, 

since the frequencies ujm+i,m and ujm,m+i are large on 
the scale of relaxational and tunneling processes. Then, 
in the case of this element can be expressed 

with the help of its dynamical equation through 



i TT 1 



(4.2) 



In the right part of this equation the terms containing 
X"", and have been dropped 

because retaining them would be against our strategy 
of going across the barrier along the shortest path to 
X™ ™. For the same reason we have also dropped the 
terms X""'^'"" and X"-'""! in the equation for X""". 
Retaining all these terms would imply taking into ac- 
count the level mixing inside the wells, which we neglect 
for small transverse fields. Now, Eq. (^^) can be it- 
erated until X"^'^^''^ is expressed through X™"^, and 
similar can be performed on X'"'™^^. Substituting their 
expressions into the equation for X™™, one arrives at the 
slow equation 



— 2 ^^"^ 



X' 



(4.3) 



where flmm' is the tunneling frequency coinciding with 
the tunneling level splitting Aem,„/ ofEq. (2.6). One can 



see now that the algorithm used here for the adiabatic 
elimination of the fast degrees of freedom in the density 
matrix equation is the dynamic counterpart of the per- 
turbative approach leading to the chain formula ( p.4| ). 
The antidiagonal matrix elements X"^"^ and X™ ™ are 
generated, in turn, by the diagonal elements X"™ and 
X"^ in ^ ^j^^ ^Y^e dynamical equations for them can be 
obtained in a similar way. The result for X™™ reads 



v-mvn ■ v'mm /-j / -r^m m v'mm\ , p 



For the matrix elements X^ 



and X' 



(4.4) 



one obtains 



equations similar to Eqs. ( [4.3| ) and ( [4.4[) . 

To formulate the resulting system of slow equations in 
a more convenient form, we introduce 



i(X'' 



X' 



These variables satisfy the system of equations 
1 



[cf. Eq. (|0)] and 



(4.6) 



where the equation firs t of Eqs. (4.7) is a con- 
sequence of Eqs. (i.6). The conservative part of 
Eqs. ( |4.7| ) describes the precession of the pseu- 
dospin (Tr,im' = {Xmm' ,Y„^m' , Zram'} in the pscudofield 
'H-mm' = {^mm' i^,^mm'}- In the absense of dissipa- 
tion, in resonance {ujmm' — 0), the pseudospin rotates in 
the y, z plane, and the difference of the level populations 
Zmm' oscillates with time. Note, however, that the Y 
and X components of the pseudospin have nothing to do 
with the actua l spin components Sy and Sx which remain 
zero, see Eq. ( 3.11 ). The only exclusion is the resonance 
between the two neighboring levels m and m + 1 near 
the top of the barrier , which is realized, e.g., for S odd 
and Hz = 0. In this case, which is actually no longer the 
tunneling case since flm,m+i oc is not suppressed by 
the anisotropy, the rotation of the pseudospin couples to 
the rotation of the real spin. 

Since the tunneling frequency ^Imm' is typically very 
small, th e co rrespondingly small detuning Umm' > ^mm' 
[see Eq. (4J)] is sufficient to suppress the resonance. On 
the other hand, a small ac field Hz{t) with a frequency 
about Qmm' giving rise to the corresponding z compo- 
nent of the pseudofield ujmm'{t) [see Eq. (4.1)] can excite 
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the tunneling resonance. The latter, however, can only 
happen under two rather severe conditions: 



H,{t) « 



(m' — m) 



(4.8) 



The former is the condition of the linear resonance, 
whereas the latter requires that the pseudospin have a 
strong preference along the x axis, in other words, that 
only the lower of the tunneling-splitted states (the even 
one) is thermally populated. The temperatures required 
by the second condition are so small that only the res- 
onance between the ground-state levels m = ±S can be 
discussed. 

The small value of the pseudofield i~lmm' in resonant 
tunneling equations (4.7) suggests an important role of 
the relaxation terms . The diagonal relaxation term R„im 
following from Eq. ( |3.13 ) has the form 



{Wm.m+lNm+l — Wm+l.mJ^m) 
+ Tfn.m-liWm,m-lN,n-l ^ W„,_ i^„,iV„) , (4.9) 

describing the exchange of particles with the levels m ± 
1. For the antidiagonal matrix elements, the relaxation 
term Rmm' in Eq. ( 3.13| ) contains X""" itself, as well as 
the matrix elements x™^^'™ These matrix elements 
do not belong, however, to the antidiagonal ones (see 
Fig. |l|); they are small slave variables that have been 
eliminated above. Dropping them leads to 



Here the terms and the analogous F^/ are the 
linewidths of the levels m and m! arising from the tran- 
sitions to the levels mil and m' ± 1 with the absorption 
or emission of an energy quantum. 

At temperatures T ^ oja = (25 — 1)-D, which is 
about 13 K for Mni2Ac, most of the particles are in the 
ground states m = ±5. The linewidths of these states 
are much smaller than that of excited ones since in Eq. 
(4.10) the emission term is absent and the absorption 
term is small as exp(— wa/^). Further lowering of the 
temperature leads to the suppression of the thermoacti- 
vation relaxation mechanism and, simultaneously, to the 
vanishing of dissipation in the ground state. Thus, the 
spin of the magnetic cluster behaves like an undamped 
two- level system (TLS). It is, however, well known (see, 
e.g., Ref. 50) that the coupling of the TLS to the bath 
strongly changes its dynamics, and one can ask where 
this coupling was lost in our calculations. The answer 



is that treating the non-Markovian relaxation term (3.S) 
we ha ve used the simplest unperturbe d t' depen dence s 
(3.12) for the spin operators of Eqs. (3.9) and (3.11), 
which do not describe the tunneling motion. This tun- 
neling motion couples, however, to a very small number 
of extremely-long-wavelength phonons, and their contri- 
bution to the relaxation terms is smaller by a factor of 



order {U-s,s/^aY e-xp{ujx/T) [see Eq. ( |3.15| )] than that 
of the regular phonon processes. Thus, the coupling of 
the tunneling mode to the bath becomes important only 
at very low temperatures. In this range serious compli- 
cations arise (see, e.g., Ref. |50| ) since the pseudospin part 
of the effective TLS Hamiltonian, TYtls = — ffi, is no 
longer large in comparison to the coupling to the bath 
and the perturbation theory breaks down. 

The equation of motion for the pseudospin, Eq. (^), 
is not closed because the relaxation term in the first line 
couples it to other levels. If we neglect t his c oupling for a 
moment, then the eigenvalues A of Eq. (4.7) determined 
as X,Y^Z (X e~^* are given by the roots of the cubic 
equation (A — F)-^A + ri-^(A — F)-|-aj^A = 0, where we have 
dropped the index mm' . This equation can be solved 
only in limiting cases. In particular, in resonance (w = 
0) the last equation of Eqs. ( ^!^ ) decouples from the 
first two ones, which describe now a damped harmonic 
oscillator with Ai,2 = (1/2)(F ± VF2 - AD?). One can 
see that the tunneling oscillations of the particle between 
the two levels become overdamped for F > 20.. In the 
small damping case, the solution of Eq. (4.7) with the 
initial condition Z{Q) = 1 has an interesting two-scale- 
relaxation form 



Zit) . exp {~^^Vt) cos(v/?^i) 

,2 



■ exp 



f72 +Cj2 



rFt 



(4.11) 



These results should not be overstated for the present 
model because in the underdamped case the neglected re- 
laxation terms in the equation for Z can be of the same 
order of magnitude as the accounted ones in the equa- 
tions for X and Y . In this case the pseudospin con cep t 
breaks down and one should use the two equations ( |4.6| ) 
instead of the first equation of Eqs. (4.7). But in the case 
of strong damping the level populations cannot deviate 
substantially from their equilibrium values because of the 
slow tunneling motion, and the different terms in the di- 
agonal relaxation terms Rmm given by Eq. (4.9) nearly 
cancel each other. Here the concept of the independent 
pseudospin is justified, and one can see that its motion 
is ind eed overdamped. Neglecting the terms X and Y in 
Eqs. (4/7), one eliminates X and Y and comes to the sim- 
ple relaxational equation for Z with A = ji^uP' ■\- F^). 

The argument in favor of the pseudospin model is that 
there can be other relaxation mechanisms, such as those 
due to spin-spin interactions, which contribute only to 
the linewidths (i.e., to the transverse relaxation rate) and 
not to the transition probabilities (i.e., to the longitudi- 
nal relaxation rate). In this typical for the magnetic res- 
onance situation the term Rm'm' — Rmm in the first equa- 
tion of Eqs. ( |4.7] ) can be neglected on the relatively short 
scale of the transverse relaxation time. In our model the 
dipole-dipole interactions could play such a role, but for 
Mni2Ac the main ef fect of such a type comes from nu- 
clear spins (see Sec. VH). 
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The possibihty of the overdamping o£lhe coherent spin 
oscillations was pointed out by Garg,LJ who considered 
resonant tunneling with the help of a phenomenologi- 
cal damped Schrodinger equation in the matrix repre- 
sentation in the unperturbed basis. Although the qual- 
itative conclusions of Garg are the same as the present 
ones, there are some discrepancies between the two ap- 
proaches in treating the relaxation. In particular, the 
eigenvalues for the two-level problem satisfy in Garg's 
approach a quadratic equation instead of the cubic or 
quartic ones in our method. Garg's solution for the 
splitted energy levels £1,2 is explicitly given by £1^2 = 
{l/2)[Ei + E2± ^{Ei- E^)^ where E, = e^- iT, 

are the damped "unperturbed" energy levels. Here the 
well-known deficiency of the damped Schrodinger equa- 
tion can be seen: The linewidths of the two levels cancel 
each other under the square root which is responsible for 
the tunneling. In the symmetric (unbiased) case this can- 
cellation is complete, and the tunneling resonance cannot 
be overdamped, in contrast to the results of the density 
matri x for malism where the linewidths are added [see 
Eq. ( 4.10 )]. This problem was avoided by Garg by con- 



sidering the resonance between the zero-width ground- 
state level in one well with an excited one in the other 
well in the low-temperature biased case, which allowed 
him to obtain plausible results. In our model in this 
case one should use Eqs. (^) with Rmm — and 
Rm'rn' ~ — 2r.,„/ iVm / , as wcll as the second and the third 
equations of Eqs. ([4.7[ ) with Ymm' = r„i', which leads 
to a quartic secular equation for A. In fact, however, 
such tunneling resonances are typically overdamped, and 
both methods give the same results. The coherent tun- 
neling oscillations should be looked for between the two 
ground-state levels whose damping is very small. For 
this situation, as well as for the description of thermally 
activated tunneling, the damped Schrodinger equation is 
inappropriate even as a qualitative tool. 

In the Arrhenius regime the rate of the process is con- 
trolled by the climbing of particles up the barrier, which 
is small in comparison to Emm' of Eqs. (4.10). In this 
case, again, one can neglect the time derivatives X and 
Y in Eqs. (4.7), which leads to the system of balance 
equations 



, 

AT- . mm' 



-{Nrn' - Nrr,) + Rmm, (4.12) 



mm' ' rjtiit 

where the rate coefficient for the transition across the 
barrier is the same as in the overdamped case and Rmm 
is given by Eq. (4.9). The form of these equations is 
quite plausible and resembling of the Fermi golden rule: 
The tun nelin g frequency fl is the transition amplitude 
[cf. Eq. dU)], whereas Ymm'/iu^m^, + E^^,) plays the 
role of the 6 function selecting the allowed resonant level 
partners. In our case of the discrete spectrum, one can- 
not set the latter to the 5 function, which causes a small 
problem: If the two levels are not exactly in resonance, 
the tunneling term preve nts e stablishing the equilibrium 
Boltzmann distribution (3.16). The corresponding devi- 



ations from the equilibrium are, however, small and they 
can be neglected, especially as we ignore all the effects 
of the level renormalization due to the transverse field . 
More important is that the tunneling term in Eq. ( 4.12 ) 



allows the establishing of the equilibrium between the 
two wells by crossing the barrier, and this process is of 
resonant character. One can speculate how the form of 
this term manifests itself in the escape rate F and what 
will be the shape of the corresponding resonances. These 
questions will be answered in the next section. 



V. ESCAPE RATE IN THE THERMALLY 
ACTIVATED REGIME 

As was said at the end of the previous section, in the 
low-tcmpcrature range T <ti U the rate of thermal acti- 
vation to the top of the barrier is much lower than that of 
the relaxation between the neighboring levels. In this sit- 
uation quasiequilibrium is promptly established in each 
of the wells, and the subsequent relaxation changes only 
the collective variables — the numbers of particles in the 
wells, N±. On this stage the problem can be solved an- 
alytically, and the solution shows that deviations from 
quasiequilibrium are localized to the narrow region near 
the top of the barrier. For the thermal activation of par- 
ticles described by the Fokker-Plank equation, this^rob- 
lem was solved in the pioneering work of Kramers .Ej The 
same method was_applied later to classical magnetic par- 
ticles by Brown.LJ For the spin system with a discrete 
spectrum the generalization was given in Ref. 1^. An- 
other method applicable in the whole temperature range, 
for small deviations from equilibrium, was suggested in 
Refs. and ^ for classical magnetic particles and in 
Ref. [42| for discrete spin systems. 

In our low-temperature case, the time derivatives in 
Eq. ( 4.12 ) can be neglected for all values of m except for 
those near the bottom of the wells, practically except for 
m = ±S. This is because the thermal activation process 
is exponentially slow and, in addition, the level popu- 
lations away from the bottoms are exponentially small. 
Now let us represent in Eq. ( 4.12| ) as 

^ 11m. , 



(5.1) 



where Nm^ is t he eq uilibrium population of the level m 
given by Eq. (3.16) and Um describes deviations from 
equilibrium. In terms of U m th e kinetic equation ( 4.12 ) 
can be with the use of Eq. ( |4.9| ) rewritten as 

— jmm' ~^ jm.m-\-l "t" Jm,m— 1: 

jmn — ^mn{,^n "^m); (^■^) 

where jmn has the meaning of the particle's current from 
the nth to the mth level, Um plays the role of a potential, 
and the conductances amn are given by 

-L mm' -^^m n^-^*^/ 



— mm 



2 cu] 



2 

mm' 



^m,m-\-l 



m,m+ 



,Wm+l.nN^\ 



(5.3) 
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where for the tunnehng process we have dropped the 
small terms violating the equilibrium Boltzmann dis- 
tribution and symmetrized the rest. One can check 
that CTm^m+i = cTm+i.m due to the symmetry of lm,m+i 
and the detailed balance condition Wm+i.mNm'^ = 

Wm^m+iN^ii- In the high-barrier hmit T < f/ the 
quantities cfm,m+i are determined mainly by the Boltz- 
mann factors and they become very small near the top 
of the barrier. On the contrary, for not too low temper- 
atures the tunneling conductances (Tmm' are extremely 
small near the bottom and increase by a giant factor [see 
Eq. (2.11)] with each step to the top of the barrier. As a 



result, (Jmm' is essential only near the top of the barrier , 
where it competes with (Tm.m+i and shunts the equivalent 
resistor circuit. 

In a broad range of m not close to either the top or 
the bottom the particle's currents jm,m-i in both wells 
are practically constant and equal to each other; let us 

denote them jm,m-i = Jh , the current from the left (— ) 

to the right well. Then one can write 



N+ = J, 



N- = -J, 



(5.4) 



for the numbers of particles in both wells. The poten- 
tial Um is also constant in the main part of the wells 
and changes near the top of the barrier where cTm.m+i 
are especially small, in accordance with the concept of 
quasiequilibrium described above. Denoting the values 
of u in the wells as m+ and u_, one can relate the differ- 
ence u+ — U- to the particle's current by the linear 

relation 

j+_ = 5-+_(u_ - (5.5) 

where is the effective barrier conductance to be de- 
termined. 

The numbers o f pa rticles in the wells, N±, calculated 
according to Eq. (5.1) are given by 



N± 



N. 



(0) 



z±/z, 



(5.6) 



where Z — Z+ + Z- is the spin partition function and Z± 
are the partition functions in each of the wells. For the 
latter it is convenient to introduce the reduced variables 



T 



S^D 
T ' 



K = 



2a 2SD' 



(5.7) 



which are equivalent to those used for the deseription of 
classical single-domain magnetic particles. e3E3 Then in 
the case of not too strong bias ft.^ ^ 1, at low tempera- 
tures the partition functions have the forms 



2 cosh^e" 



1 



-2a/S ■ 



1 



-2a/S ■ 



(5.8) 



Combining now Eqs. (x4), (5_^), and ( |5.6| ) one comes to 
the rate equations 



N+ 



N± 



N- 



(0) 



N. 



(0) 



(5.9) 



For the average spin polarization 

m, = {S,)^S{N+-N^) 
the latter result in 



^(0)^(0)^ 



(5.10) 



(5.11) 



and (p|), 7Vj°)iVi"' 



rhz — ~T{mz — ml 

where, according to Eqs. 
(4cosh2 0"^ 

Finding the effective barrier conductance deter- 
mined by Eq. ( pTs] ) is the easiest task in the case without 
a transverse field where amm' = 0. Here the elementary 



resistances cr,m „j_|_i of Eq. (5.3) add with the result 



s-i 



"H — ^ 2^ '-'m,m+l- 



(5.12) 



m=-S 

For the thcrmoactivation rate F this yields 

-, -1 



F ^ 



4 cosh^ ^ 



E 



exp{em/T) 



(5.13) 



One can see that the main contribution to this expression 
comes from the top region, so that F cx exp[— a(l — h^)'^] 
and t he ex act limits of summation i n Eq s. ( |3.12 ) and 
Eq. ( 5.13| ) are irrelevant. Formula ( [5.13D ri^ the micro- 
scopic generalization of the Brown's rcsulto on systems 
with a discrete spectrum. For S — ]_a similar result was 
obtained in early work by Orbachjla and for a general 
spin generalizations were given in Refs. ^ and 42 in the 
unbiased and biased cases, correspondingly. In R cf. |4^ 



different limiting forms of the prefactor in Eq. (5.13) 



were analyzed. The most striking of its features is its de- 
pendence on the bias field Hz with a strong decrease in 
the region where two levels at the top of the barrier come 
into resonance. The latter is due to the frequency depen- 
dence ( 3.15| ) of the one-phonon transition rate between 
these levels, Wm+i,m- 

In the case of a nonzero transverse field the barrier 
conductance can be calculated by a well-known re- 
currence procedure starting from the top of the barrier. 
Introducing amm' as the total conductance due to the 
part of the barrier between the "points" m and m' (see 
Fig. nf) one obtains 



1 



m+l,m' 



m.m-\-l 



(5.14) 



with a proper initial condition at the unperturbed top 
of the barrier , m„iax ~ Hz/ {2D). If the spin is large 
and the transverse field H^ is not too small, the level 
pair TOf), m'f, corresponding to the actual renormalized top 
of the bar rier is situated many "steps" below mmax [see 
Eq. (2^)]. In this case the starting point TOrnax be- 



comes unimportant, and the recurrence algorithm (5.14) 
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generates a continued fraction. In the Arrhenius regime 

T <^ U, the quantity CTmm' rapidly converges to 

down from the renormahzed top of the b arrier rnh,m[,. 



The role of different terms in Eq. ( 5.14 ) can be made 
clear if one considers the ratio 



2 

mm' 



^m+l,m 



^m+l,m 



r2 , 

mm 



li^mmH ^ li^m+ljml. 



(5.15) 



corresponding to the nonresonant and resonant situa- 
tions. If this ratio is of order unity for some pair to;,, to'^, 
one can consider all the tunneling conductances LUmm' 
above thi s lev el as infinite and below this level as zero 
[see Eq. (2.11)]. In the resonant situation, one also can 



speak about conducting and blocked resonances. Since at 
the level to& -I- 1, — 1 the circuit is completely shunted, 
one concludes that renormahzed by the transverse field 
the top of the barrier is localized at m — mi,, with an 
uncertainty of one level. In the nonresonant situation for 
1 <C |m| ^ S* this leads to the previously obtained clas- 
sical result of Eq. (2_^). At resonance, for Hz = 0, the 
corresponding value of is determined by the equation 



2S^K = 



2D\mb\ 



(5.16) 



Since the level linewidths are small, T^m' ^ D, this 
value of TOb is greater than that off resonance, which 
thus leads to the resonant dips in the effective barrier 
height. Note, however, that the magnitude of these dips 
i s str ongly reduced by the exponent l/(2|m6|) in Eq. 
(5.16), so that they become small in systems of large 
s pin. The shape of resonances in the escape rate F of Eq. 



(5.11) can be visualized, if one considers resonant transi- 
tions between only one pair of levels to^jtoJ^. Neglecting 
transitions above this level, one writes 



1 



(5.17) 



the resonance broadens away. But there are tunneling 
resonances between lower pairs of levels for which the 
same formula ( ^.18 ) can be written. The width of these 
peaks ilm,m' is much smaller, but their height at reso- 
nance '-^ Tm^m'Nm'' incrcascs with the level depth as the 

Arrhenius factor iVm^ ''-^ exp(— Em/T) and is maximal for 
the deepest unblocked pair of resonant levels. In fact, in 
the low-damping case t he li ne shape of T described by 
the continued fraction ( ^.14| ) consists of many peaks of 
stepwise decreasing width ~ ^m.m' mounting on top of 
each other and forming a self-similar structure. 

An illustration of the behavior of the escape rate F in 
the Arrhenius regime based on numerical calculations of 
the barrier conductance will be given in Sec. VI] . 



In the next section we briefly discuss the range of lower 
temperatures where a "more quantum" behavior of F is 
to be expected. 



VI. TUNNELING VERSUS THERMAL 
ACTIVATION 



In the Arrhenius regime above, the product fl^^iN, 



(0) 



in the tunneling conductance <Tmm' of Eq. (5.3) increases 
unlimitedly up the barrier and amm' shunts the effective 
circuit at some level mb determining the renormahzed 
position of the top of the barrier. This mechanism is of 
resonant character, but the temperature dependence of 
the escape rate remains classical. With lowering tem- 
perature the question arises, of which group of levels the 
tunneling conductance <7mm' has a maximum. The anal- 
ysis of the function /(m) = exp(— Em/T) shows 
that there are two more regimes in addition to the Ar- 
rhenius one — ground-state tunneling and thermally as- 
sisted tunneling. The temperature of the crossover be- 
tween these two regimes, Tqo, is determined from the 
condition f{—S) = f{—S + 1); i.e., the rate of tunnel- 
ing from the first and other excited states falls below the 
ground-state tunnelin g rat e. The value of Tqo calculated 
with the help of Eq. (2.11) has the form 



where cr__ „^ is the conductance between the bottom of 
the left well and the point mh, etc. This expression can 



be rewritten with the use of Eq. (5.3), and for the escape 
rate F one obtains 



mb ,mjj 



AO) 

•Imt 



F2 , 



An'' 



(5.18) 



where A 



F >N^°Ua 

-L mf, ,m'^-^^ mi, yJ - 



-1 



-(T-^l+)- FromEqs. (|5 



I) 

and (4.10) it follows that A ~ 1, if the resonant transi- 
tions through the lower- lying pairs of levels are neglected. 
Thus, contrary to what could be naively expected, the 
linewidth of the resonance in the escape rate F is in- 
sensitive to the level linewidth T^af,,m' which is smaller 
than the tunneling frequency fim^^mj^ ^or conducting res- 
onances. This frequency grows rapidly with the trans- 
verse field. When it reaches the level spacing |ijj„i+i,m|, 



Tn, 



SD 



He'S/hl) 



2SD' 



(6.1) 



In theories of tunneling using continuous level models the 
quantity Tqo does not appear. For models with discrete 
levels one should keep in mind that the linewidth of the 
ground states is much smaller than that of excited ones, 
and this should make the analysis more subtle, but we 
will not further pursue this topic here. 

For T > Too tunneling goes through the group of levels 
between the bottom and the top for which /(m) has a 
maximum; if the position of this group does not coincide 
with the top of the barrier to = to^, this regime is called 
thermally assisted tunneling. There are different scenar- 
ios for the temperature dependence of this group of levels, 
TO ^ totat- It can shift continuously from the bottom to 
the top with a crossover to the Arrhenius regime at some 
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FIG. 3. Temperature dependence of the group of lev- 
els, rriTAT, making the dominant contribution into the ther- 
mally assisted tunneling, determined from the maximum of 
/(m) = i^^^, exp(— Em/T) in the unbiased case. 



temperature Tq. The other type of behavior is reahzed if 
the function /(m) has two maxima, say, at the top and 
near the bottom of the barrier. In this case there are two 
competing channels of relaxation which go from one into 
the other at the crossover temperature Tq. Both of these 
scenarios were studied for the models with continuous 
spectra, and the analogy with the second- and first-order 
phase transitions was pointed outE3 

For the uniaxial spin model both types of thermally 
assisted tunneling can be realized, and the situation can 
be controlled by the transverse field. In particular, for 
the second-order transition the cross over temperature Tq 
obtained with the help of Eq. ( ^.11 ) is given by 



(2) 



8 2 

(2) 



(6.2) 



For low transverse fields Tq ' becomes too small, and 
the first-order transition to the regime of thermally as- 
sisted tunneling occurs when the temperature is lowered 

(2) 

before Tq ' is reached. Details of the analysis will be 
presented elsewhere; here we illustrate the temperature 
dependence of m ~ ?titat in Fig'- It can be seen 
that the higher values of hx favor the second-order tran- 
sition: The curve mTATC?") goes "continuously" through 
each value of m and merges at Tq with the horizontal 
line TO = TOb characterizing the Arrhenius regime. On 
the contrary, in lower fields hx large jumps of totat at 
To can be seen. For smaller spins the low-temperature 
tail of the curve totat(T) becomes shorter. The value of 
Too is in all cases well described by formula (|6J|). 

In the thermally assisted tunneling regime, the ratio of 
the tunneling and intrawell conductances, Eq. ( ^.15 ) is 
a very small number in the relevant region m ~ totat- 
Thus the slow tunneling process controls the escape rate 
r, and the distribution of particles in the wells does not 



deviate from quasiequilibrium. In this case F is simply 
given by 



F = 



^(0)^(0)' 



(6.3) 



va— — S 



i.e., it is the tunneling probability weighed with the 
Boltz ma nn factor [see Eq. (5.3)]. Expressions similar to 
Eq. (6.3) were taken as a starting point in many investi- 
gations of the escape rate of particles from a metastable 
well at nonzero temperatures (see, e.g., Ref. 57). An 
efficient method of treating this problem for continuous 
spectra, including the-dissipative case, is based on the 
instanton techniqucuS For our spin model, however, 
the spectrum cannot be made continuous by a reason- 
able variation of some physical parameter; the tunneling 
frequency changes abruptly from one level to another, 
and this situation persists in the limit S —^ oo (see the 
end of Sec. |l|). This situation seems to be pertinent 
not only to spin systems, which can be, in fact, mapped 
onto the particles ,l3'c3 but for double- well models in gen- 
eral. Resonant tunneling between the discrete levels in 
a low-damped SQUID was observed recently in Ref. |5^. 
The numerically calculaJxid tunneling level splittings for 
the SQUID HamiltonianE3 also change abruptly from one 
level pair to another. 

The advantage of our more general approach to find- 
ing the barrier conductance cr^ based on the recurrence 

relations (5.14) in comparison to the simplified formula 
( |6.3D is its ability to handle the case of very small cou- 
pling to the bath. In this case the relaxation rates for 
the exchange between the neighboring levels Um,m+i of 
Eq. (5.2) become very small, as well as the tunneling 
conductances amm' off resonance, and so does the result- 
ing escape rate F. If one sets the system on resonance 
to increase tunneling, then the system does not co me to 
quasiequilibrium in each of the wells and formula (6.3) 
breaks down. 



VII. NUMERICAL RESULTS FOR THE ESCAPE 
RATE; ROLE OF NUCLEAR SPINS AND THE 
AXIS MISALIGNMENT 

In this section we present the results of numerical sim- 
ulations for the escape rate F obtained with the methods 
of the previous section in the Arrhenius regime. The re- 
gion below the crossover temperature Tq is not further 
considered in this paper. For systems of moderate spin 
the range of thermally assisted tunneling is rather nar- 
row, and at temperatures T < Tqq in the unbiased case 
tunneling should go between the ground states. Since 
the linewidths of the ground-state levels are exponen- 
tially small at such temperatures, even a small detuning 
is sufficient to suppress the resonance. In this case we 
face a strongly nonresonant situation, and our theoret- 
ical methods of Sec. IV should be modified. Even in 



the Arrhenius regime, there is a problem with nonres- 
onant processes — the escape rate calculated with the 
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FIG. 4. Bias-field dependence of the resonant tunneling es- 
cape rate of the uniaxial spin model. 



help of Eq. ( 5.14 ) shows discontinuities at values of the 
bias field, at which we switch from one resonant partner 
level to another in the calculation routine. In fact, the 
resonance of each level with several partners in the other 
well should be considered, but a rigorous treatment of 
this problem would lead to serious complications. For 
this reason we simply extend the applicability of the ki- 
netic equation ( 1.12| ) by considering, for each level to, 
the tunneling resonances with the two partners to' and 
to' — 1 satisfying e^' < £m < £m'-i- In this symmet- 
ric approach the switching between partners occurs in 
resonance and no discontinuity in V appears. The calcu- 
lations in this case can be performed wit h the help of the 
modification of the recurrence relation ( 5.14 ). The reso- 
nance oLone level with all other partners was considered 
by GargU using the damped Schrodinger equation. 

Treating the re lax ation terms we replace lm,m+i by 
lm,m+i [see Eq. (3.8)], which amounts to dropping the 
operator Sz in Eq. (3.4). Then we fit the strength of the 



spin-phonon coupling to the measured for Mni2Ac value 
of the prefactor Fq = 5 x lO^s""'^ in the escape rate F = 
Fq exp(— t/eff/r). Since the experimental temperatures 
about several kelvin exceed the level spacing near the 
top of the barier, uim+i,m ~ 2Dm ~ 1 K, the prefactor 
depends linearly on temperature; see the first line of Eq. 
( |3.15 ). This dependence is, however, difficult to see in 
the limited temperature interval. 

The results for the escape rate as a function of the bias 
field Hz are represented on Fig. ^ for different values of 
the transverse field Hx- One can see the superpositions 
of broad and narrow peaks at the resonant values of the 
bias field Hzk = Dk, which correspond to the tunneling 
via the shallower and deeper resonant levels, respectively. 
The width of peaks alternates as a function of the reso- 
nance number k, since the tunneling transitions between 
different pairs of levels appear in even or odd orders of the 
perturbation theory in Hx/D; see Sec. ^ In particular. 



in the unbiased case fc = 0, the tunneling between the 
level pair —1, 1 appears in second order, il-i^i oc H^/D. 
As a result, a very narrow peak in F emerges at Hz — 
for hx = 10"'*. This peak broadens with the increase 
of Hx, and at hx = 5 x 10'^ a new narrow peak corre- 
sponding to the resonance —2, 2 with « H!j:/D^ is 
seen. A similar picture holds for k — 2 and other even 
resonances. 

For the odd resonances, as fc 1, the escape rate 
in zero transverse field becomes small due to the fre- 
quency dependence of one-phonon processes discussed 
above. The same result was obtained for the tunneling 
assisted oae-phonon processes between the deep levels in 
the wells.El In Fig. we have included a small frequency- 
independent contribution from the Raman scattering 
processes to obtain a nonzero value of the escape rate.cHI 
This feature is, however, completely suppressed already 
in very small transverse fields because of the opening of 
a new transition channel: the tunneling between the top- 
most resonant pair —1,0 that appears in the first order, 
ri_i.o cx: Hx- The latter is, in fact, a kind of a free pre- 
cession around H^; rather than tunneling. The rate of 
this precession competes with the small relaxation rate 
[see the second line of Eq. ( 5.15| )]; that is, the purely 
dynamical transition between the levels —1,0 competes 
with the dissipative one. As a result, the dip in F yields 
to the massive peak already for = 10^* for the damp- 
ing parameters appropriate for Mni2Ac. 

For higher values of the transverse field h.^ , the behav- 
ior of the even and odd resonance peaks is the same. As 
hx is growing, the condition ^mm'{hx) > ^mm' of Eq. 
(^T5|) for a given pair of levels to, to' is satisfied at a 
certain value of the transverse field hxt- At that value 
the TO, to' resonance becomes unblocked. This results in 
a new narrow peak, of width about ^mm' (hx) , which ap- 
pears on the top of the to + 1 , to' — 1 resonant peak (see 
Fig ^) . This situation is quite universal in the sense that 
the TO-t-1, to' — 1 resonant peak can itself be a narrow peak 
on the top of the to -t- 2, to' — 2 resonant peak. In general, 
each resonance consists of a few peaks mounting on top 
of each other. The width of two consequent peaks within 
one resonance differs by a factor about ~ 55, in accor- 
dance with Eq. ( 2.11 ). The magnification of the Hz inter- 
val around the resonant values Hzk shows the self-similar 
multitower structure of the resonance. In that structure 
the total number of peaks depends on the strength of the 
dissipation, while their height is determined by tempera- 
ture. The lower the damping, the greater is the number 
of the peaks. The lower the temperature, the greater is 
the difference in the height of the peaks mounting on top 
of each other. 

The dependences of F on the transverse field Hx for the 
resonant and slightly off-resonance values of the bias field 
are shown in Fig. |[ The steps on the resonant Hx depen- 
dences of F correspond to the valu es of Hx at which the 
value of TOf, determined from Eq. (5.16) takes an integer 
value (or a half-integer value for systems of half-integer 
spin S) . For these values of Hx a resonant shunting of the 
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FIG. 5. Transverse-field dependence of the resonant tun- 
neling escape rate of the uniaxial spin model. 



barrier at the next deeper level occurs. The flat regions 
correspond to the situation when one pair of resonant lev- 
els is already completely shunted and the following (the 
lower) one is yet completely unshunted. The step val- 
ues of Hx are sensitive to the sum of the level linewidths 



given by Eq. ( 4.10| ), and thus such experiments 



are conceivable as a kind of spectroscopy measuring the 
relaxation characteristics of separate levels. The id plot 
of T{HxtHz) summarizing the features of resonant tun- 
neling process discussed above is presented in Fig. |^. 

Apart from resonant tunneling, the overall shape of 
T{Hxj Hz) follows approximately the Arrhenius law with 
the classical barrier height U{Hx, Hz). This can be seen 
especially clear for systems with large spin, frequency- 
independent relaxation rates and low temperatures. The 
last condition is needed to reduce the relative role of the 
field dependence of the prefactor Fq = To{Hx) in the clas- 
sical expression for F, which is not yet well established 
(see Refs. |^ and ||). The comparison of our calculation 
with the classical result accounting only for the depen- 
dence 



UiHx,Hz) = DS''il-hzf 



1 - 2h. 



(7.1) 



for /ij; <C 1 is presented in Fig. |^. The rather good 
accordance between the classical and quantum results il- 
lustrates the conjectures of Sec. || in a more general 
biased case. 

The resonant tunneling curves obtained alMyejidfl-apt. 
fully explain the experimental observationscjOEaO'EHl 
showing that all peaks have approximately the same 
form. The latter can be the consequence of the averag- 
ing effect due to the misalignment of the particle's axes in 
not perfectly oriented polycrystalline samples. A similar 
effect can be caused in Mni2Ac by nuclear spins whose 
fluctuating transverse components can, in addition, in- 
duce tunneling even in the absense of an externally ap- 
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FIG. 6. Dependence of the escape rate of the uniaxial spin 
model on the field components Hx and Hz (parameters are 
appropriate for Mni2Ac). 



plied field Hx- The corresponding adjustments of our 
method will be made below. 

In a Mni2Ac molecule each of 12 Mn atoms interacts 
with its own nuclear spin 1^, / — 5/2, via the hyperfine 
(HF) interaction. For the total cluster spin this interac- 
tion can be approximately written as 



-AoffSItot, 



It 



(7.2) 



In fact, the hyperfine interactions are somewhat differ- 
ent for different Mn atoms, and their extensive discus- 
sion can be found in Ref. If all the nuclear spins 
are aligned in the same direction, the energy of the 
HF interaction i?HF,max — l2ISAcs — 0.6 K is com- 
parable to the level spacing near the top of the bar- 
rier Iwm+i.ml ^ 2Dm ^ 1 K and is much greater than 
the dipole-dipole energy ~ 0.06 K. The effective 
HF field produced by the nuclei on the cluster spin is 

in this case about i?HF,max — £^HF,max/(5'AiBS') ~ 0.05 

T. If this HF field is perpendicular to the easy axis z, 
the corresponding dimcnsionlcss transverse field hx^uF — 
gfiBHiip^jnax/{2SD) ~ 3.7 x 10^'^ should result in stron, 
resonant (as well as nonresonant) tunneling; see Fig, 
On the other hand, the role of the z component of the 
HF field in resonant tunneling is determined by the di- 
mensionless parameter gpLBHBF,-ma:^/ D ~ 0.075. This 
shows that the narrow resonance lines in Fig. ^ should 
be averaged away by the fluctuating z component of the 
hyperflne field; i.e., the hyperfine interaction suppresses 
the coherence^ This second effect was discussed by sev- 
eral authors Jl3 here we will take into account both effects 
of nuclear spins with the help of simplified qualitative ar- 
guments. 

The subtlety of the hyperfine interaction is that it 
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FIG. 7. Longitudinal field dependence of the escape rate 
of a large-spin uniaxial model in the fixed transverse field. 
The curve 1 was obtain ed f rom the pure thermoactivation 
curve fea; = using Eq. (7.1) for the barrier lowering in the 
transverse field = 0.1. 



conserves the total projection 5*2 + ^iliz, and, strictly 
speaking, the coupled equations of motion for the tun- 
neling cluster spin and rotating nuclear spins should be 
solved. In the Arrhenius regime, however, tunneling oc- 
curs near the top of the barrier where it is rather fast 
— it ranges from f2mm' ~ |t^m+i,r?i| off resonance to 
f^rnm' ~ ^mm' at rcsonancc. This is much faster than 
the nuclear relaxation rate which is due to the fluctu- 
ating magnetic fields and is determined by the small nu- 
clear magnetic moment. Further, tunneling of the cluster 
spin near the top of the barrier leads to a relatively small 
change of its z-projection: ASz ~ 2m{, <C 25. This is not 
a large part of the whole integral of motion Sz + J2i ^iz- 
Indeed, for the randomly oriented nuclear spins the sec- 
ond term of this sum is on average of order \/T2/ ~ 8.7, 
and thus tunneling of the cluster spin can be compen- 
sated by the corresponding rotation of the nuclear spins. 
(On the contrary, for tunneling from the ground state 
at T < To the z projection change is AS'z — 2S = 20, 
and this process cannot go via the interaction with the 
nuclear spins — it is blocked by the conservation law.) 
Thus, in the Arrhenius regime one can qualitatively con- 
sider nuclear spins as frozen — they do not change their 
state as a result of the tunneling of the cluster spin. The 
distribution function of the HF field on the cluster spin 
can be easily found. As the energy of the interaction 
of one nuclear spin with the cluster spin ISA^s ~ 0.05 
K is much smaller than temperature, one can use the 
infinite-temperature distribution function for the individ- 
ual nuclear spins. Then, for a large number of nuclear 
spins, = 12 ^ 1, the quantum-statistical averages of 
the total nuclear spin Itot in Eq. (7.2) are given by the 
Gaussian distribution function 
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FIG. 8. Resonant tunneling escape rate in Mni2Ac due to 
nuclear spins and the crystallite misalignments. 



F(Itot) 



1 



(27rfT/)3/2 



exp 



^tot 



2a,) 



(7.3) 



where the dispersion ct/ — {N/3)I{I + 1) can be checked 
calculating the average (/tot,z) directly and from Eq. 
(7.3) and comparing the results. Now, all the previously 
obtained expressions for the escape rate F, as well as such 
quantities as the time dependence of magnetization and 
dynamic susceptibility, should be averaged with the dis- 
tribution function F. In the absense of an externally ap- 
plied field Hx, the averaging of each quantity A{Hx; Hz) 
is done explicitly as 



A{Hx,Hz) ^ jdx 2x6'''^ jdz l-^A{xHuF,Hz + zHiiF), 





Hup — i^HF.max- 



(2^7) 
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(7.4) 



The results of this averaging for the escape rate F are 
presented in Fig. ||. The role of nuclear spins in induc- 
ing resonant tunneling and suppressing the narrow res- 
onance lines is clearly seen. In addition, we have taken 
into account small fluctuations of the directions of the 
anisotropy axes of Mni2Ac molecules in polycrystalline 
samples with the dispersion of only = 3°. These mis- 
alignments also produce a static fluctuating components 
of the transverse field, and their role becomes progres- 
sively more important with the increase of the bias field 
Hz- One can see that when all these effects are taken 
into account, all the resonant tunneling peaks become 
approximately of the same form, as observed in experi- 
ments. 



VIII. DISCUSSION 

We have presented the theory of thermally activated 
resonant spin tunneling. The bulk of the theory applies 
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to any molecular magnet, while particular numerical il- 
lustrations were made for Mni2Ac. Quantization of spin 
levels, which is the key to explaining experimental re- 
sults, has dictated our choice of theoretical apparatus. 
Rather than employing instanton methods, suitable for 
models with continuous spectra, we have used the density 
matrix description of the spin interacting with thermal 
bath. 

In continuous models three regimes for the escape rate 
r are usually studied. At high temperatures quantum- 
mechanical effects are not important, and the escape over 
the barrier is due to pure thermal activation described by 
the Arrhenius law. In the limit of zero temperature only 
tunneling out of the ground state is important. There 
is also an intermediate regime which combines thermal 
activation to excited levels with tunneling across the bar- 
rier, which is called thermally assisted tunneling. In that 
regime the position of the narrow group of levels which 
dominate the escape rate depends on temperature, mov- 
ing continuously from the ground state at T = to the 
top of the barrier at the temperature called the crossover 
temperature. This situation describes a conventional, 
smooth, second-order transition from quantum tunneling 
to thermal activation.Ell In principle that transition can 
also be first order, which would correspond to the sharp 
crossover frompquantum tunneling to thermal Arrhenius- 
type behavior E3 We have demonstrated that this is ex- 
actly what happens for a spin system in low transverse 
field. Correspondingly, the experimental study of the es- 
cape rate should find the evolution from sharp to smooth 
crossover between thermally assisted tunneling and the 
Arrhenius regime on the transverse field. 

In systems of moderate spin, such as Mni2Ac, ther- 
mally assisted tunneling occurs in a rather narrow tem- 
perature range. In experiments the Arrhenius law that 
occurs in a wider temperature range has been observed. 
Despite the purely classical temperature dependence of 
the relaxation in the Arrhenius regime, the field depen- 
dence of r shows quantum effects due to the discrete 
nature of spin ignored in continuous models. Contrary 
to these models, which start with a given barrier, a well- 
defined barrier does not exist for a mesoscopic spin; its 
effective value depends on the bias field in a non- 
monotonic manner. The observed minima of the effective 
barrier are due to the crossing of the spin levels, which 
results in resonant tunneling between the wells. This is 
different from a classical spin system where the barrier 
monotonically decreases with increasing ■ This regime 
can be called thermally activated tunneling, as different 
from the regime of thermally assisted tunneling. The dif- 
ference between the two regimes is that in the first regime 
tunneling always occurs at the top of the barrier, while in 
the second regime it occurs from excited levels between 
the bottom and the top of the barrier. 

The theory predicts that each resonance in the escape 
rate T has a multitower structure with peaks of decreas- 
ing width mounting on top of each other. This effect is 
due to resonant spin tunneling between different match- 



ing levels. All peaks are centered at the same field, if 
the corresponding pair of levels match at the same value 
of the bias field. Note that this assumpion relies on 
the simple form of the Hamiltonian used in our calcu- 
lations. Additional terms of different symmetry would 
violate this assumtion. If these terms are small, as they 
are in Mni2Ac, the resonances on will not be exactly 
equidistant and the centers of peaks towering in each res- 
onance must be slightly displaced with respect to each 
other. The number of peaks in each resonance increases 
with decreasing dissipation. 

Depending on the number k of the resonance [see Eq. 
(p. 2D], the leading contribution to the rate appears in 
even or odd orders of the perturbation theory on the 
transverse field . This results in the alternation of the 
shape of resonances on Hz ■ Another effect predicted by 
the theory is the stepwise dependence of the rate on the 
transverse field when the longitudinal field is tuned to 
the resonance. 

The origin of the terms in the Hamiltonian responsi- 
ble for tunneling is different for different molecular mag- 
nets. The absense of any selection rules for resonances in 
Mni2Ac unambiguously points to transverse fields caus- 
ing the transitions. These fields originate from the hy- 
perfine and, to a smaller degree, from the dipole-dipole 
interactions. In Feg the hyperfine interactions are negli- 
gible, and the transitions are presumably caused by the 
transverse anisotropy. 

Our theory for Mni2Ac can pretend to the quantitative 
description of the magnetic relaxation in this system, as 
it takes into account all major contributions to the ef- 
fect. However, observation of more subtle effects, such as 
the multitower structure of resonances, the alternating 
shape, stepwise dependence on the transverse field, etc., 
is less likely in Mni2Ac. This is because of the smearing 
of these effects by strong fluctuations of the hyperfine 
field. Feg (see, e.g., Ref. ^2|) seems to be a better candi- 
date for observing these effects. 
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